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Abstract 

The time dependent Schrodinger equation for two identical and charged pi- 
ons is solved using wavepacket states. It is shown that the expected Coulomb 
distortion in the momentum correlation function is obliterated by the dispersion 
of the localized states, and therefore becomes unobservable. 
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1 Introduction 



The Bose-Einstein correlations of identical charged pions are believed to be distorted 
by the Coulomb interaction between both pions. Since the end of the seventies, several 
efforts where made to investigate these distortions. A reliable information about the 
size of the pion source can only be extracted from the correlation signal when the 
properties of the Coulomb distortions are sufficiently understood. The size of the 
Coulomb effect, and how it should be calculated, was always a source of intensive 
discussions. In the early theoretical approaches, the square of the relative Coulomb 
wave function for a pointlike source, the Gamow factor, was employed to obtain the 
Coulomb correction of the correlation function [^. Later it was realized that this 
ansatz overestimated the Coulomb suppression, and since then different models where 
developed which successively led to smaller and smaller distortions. First, the finite 
source size was included |@] , later the infiuence of the exchange terms , the infiuence 
of long-lived resonances [Q, of Coulomb screening due to other charged particles 
and the finite momentum resolution of the detector [^]. 

Already in 1986, however, it became clear that a correct treatment of pions in the 
fireball must include their localization in terms of quantum mechanical wavepackets, 
the spatial widths of which are connected to the mean free path of the pions which has 
a value of a few fm [0]. Recently, a consistent wavepacket description was presented 
by the authors and applied to different numerical codes [§. A somewhat different 
ansatz, which may be suitable especially for very high multiplicities, is developed in 

0. It was demonstrated, how the localization of the pions introduces a zero point 
energy, which consequently leads to a momentum uncertainty in the pion ensemble, 

1. e. to dispersion. In order to estimate the possible influence of the dispersion onto 
the Coulomb effects, we compare the size of the Coulomb potential with the zero point 
energy using a simple picture: Imagine two pions which are separated by the distance 
Tq. Both states shall not overlap significantly, so that we require them to be localized in 
Gaussian wavepackets with the spatial extension cxo = ^-0/2. We evaluate the Coulomb 
energy -Ecoui = e^/ro and the zero point (quantum-) energy = 3Ti^ / {8ma'^) for a 
distance of Tq = 4 fm and obtain -Ecoui ~ 0.4 MeV and Eq ^ 26 MeV. Therefore, 
the quantum energy E^ is nearly two orders of magnitude larger than E'coui- Only for 
distances > 290 fm exceeds the Coulomb energy the localization energy and the 
picture of two classical particles with Coulomb repulsion becomes valid. For smaller 
distances the dispersion introduces a strong random motion into the ensemble which 
obliterates the directed motion caused by the Coulomb repulsion. This suggests that 
a correct quantum mechanical treatment of the Coulomb interaction is necessary and 
can not be substituted by any semiclassical model. 
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2 Matrix representation of the Schrodinger equa- 
tion 



The Schrodinger equation for two identical charged bosons, 



with the Hamiltonian 



H = -|^(Ai + A2) + I I, ^f + V (2) 

2m ri — r2 



and the symmetrized 2- particle state 

|^(ri,r2,t)) = ATs (|xi(ri,t)) |x2(r2,t)) + |X2(ri,t)) |xi(r2,t))) , (3) 

has to be solved. Since the Hamiltonian does not depend on time, the solution could 
be found on principle by time independent methods. The problem, however, are the 
initial conditions of the system: At t = each particle is localized within a certain 
volume (wavepackets in the fireball). There exists neither an analytical solution for this 
problem nor a way how to transform these initial conditions into appropriate border 
conditions for a time independent treatment. We therefore have to solve the time 
dependent problem. To start with, we expand the state (|^) into a set of basis functions 

|^(ri, r2, t)) = Ck{t) |^fc(ri, r2, t)) (4) 

with the c-valued coefficients 

Ckit) = (^fe(ri,r2,t)|^(ri,r2,t)) . (5) 

Now the Schrodinger equation can be written in matrix form 

ihc = S~^ (H-D) c (6) 



with the time dependent coefficient vector c, the Hamilton matrix Hij = (\E'j|T+V^|\E'j), 
the overlap matrix Sij = (\E'j|\E'j) and the time evolution matrix Dij = (^jjifi^t |\E'j). 
This representation of the Schrodinger equation is equivalent to Eq. (|l|) if the basis 
provides a complete set of states in the 2- particle Hilbert space. In any numerical 
calculations, of course, only a limited number of basis functions can be used so that an 
adequate choice of the basis is essential in order to reach the required accuracy and, at 
the same time, to keep the numerics within affordable limits. 

3 Choice of the basis set 

For the single particle states we use the following time dependent wave functions: 



X,>(r, t) = exp (P,(t) ■ r - 9r^^,it) ) - " ^^^^^^^ J (7) 



' (t. „ n u^\ (r-RfcW) 



2^ 
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with the phase 



ft I Tit \ 

Of^ j^it) = El'^'ir) dr + h + Ijy + + 3/2) arctan ( 1 (8) 

and x''f^'° = — Rfc^)'^"", and similar for the y- and z coordinates. The capital co- 
ordinates (Rfc(t), Pfc(t)) denote the phase space coordinates of the fc'th wave packet 
centre, whereas r is the free variable of the state in configuration space representation. 
The time dependent width is given by s{t) = cTq + iTitj {2mao) and the initial width is 
o"o. With the index vector Ij = (0, 0, 0) we get the usual Gaussian wave packet which 
was already used in ^j. Following the notation used in the harmonic oscillator prob- 
lem, this function may be called "s- type" function, and consequently the index vector 
Ij = (1, 0, 0) defines a "p- type" function and so on. We now build symmetrized 2- 
particle states from the single particle states, because then the basis functions already 
fulfill the correct permutation symmetry. An additional symmetry requirement is that 
the squared amphtude must be symmetric with respect to the centre of mass of the 
system, which requires the index vectors of both single particle states to be equal. We 
therefore define the basis in terms of the homogeneous states 

~ xra(ri, t) Xi^2i^2, t) + xr^2(ri> ^) xra(r2, t) ■ (9) 

We use a set of 2- particle basis functions, which contains 1 s- function, 3 p- functions, 
6 d- functions and 10 f- functions. The use of even larger basis sets is prohibited by the 
exponentially increasing numerical effort. For example, the evaluation of a complete 
trajectory (which typically requires 500 timesteps) needs 1 minute for a sp- basis, 1/2 
hour for a spd- basis and 1/2 day for a spdf- basisQ. 

4 The solution of Schrodinger's equation 

In order to solve the Schrodinger equation (^), the matrices V , T, D and S have to be 
evaluated. It is possible to obtain analytical expressions for the kinetic energy- and time 
evolution matrix elements, whereas the potential matrix elements, representing a 6- fold 
integral, can be reduced to an one dimensional integration and tabulated. Details of 
the very lengthy calculations may not be given here. We refer to the literature in 
quantum chemistry^. 

At the beginning of each event, the system has to be initialized. As source function 
we used the following simple parametrization of the phase space: 

''W-^^-p(^) /(P) = exp (5^) . (10) 

Values obtained on a DEC 250^/266 workstation. 

^For example, in Chapter 6, the calculations are done for the time independent Gaussians. 
For the evaluation of our time dependent basis set, a detailed script (13 pages) can be obtained from 
the authors. 



3 



i.e. the source is a static Gaussian with rms- radius y3/2Rs and temperature T. 
This choice is motivated by the fact that an analytical expression for the 2- particle 
correlation function is available for noninteracting wavepackets 

oo / 2h2 \3/2 

C2(Pi,P2) = E(-l)' l^Jf^) (e^'=(P^'P^) + e^'=(P-P^)) (11) 



k=0 



with 



/fc(Pi,P2) = -(pi-P2)'T7r^ , (12) 




and 



5fc Pi,P2 = - P1-P2 09 6^ ^ + of,2p2A ^ ' ^3 

SZa'^akCk 2h Ribkdk J 



1 1 , 1 

ak = TT^ + 7:-^ + -r^ , Ofc = + 



2(t2 2mT Aal ' i?^ 4^2 
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Cfc = at- , dk = hk- , (14) 

IGa^afc 16cr46fc 

(jp = %/ (2cTo) is the width of the wavepacket in momentum space. The state was initial- 
ized to be a pure s- state, i. e. with the coefficient vector c(t = 0) = (1, 0, . . . , 0) in Eq. 
(^), and the initial width of the Gaussian was CTq = 1.8 fm. The time propagation was 
performed using a second order differential scheme [|1T|. For each timestep, the force 



between both wavepackets was obtained numerically using P2 = — V2i(^E'|V^(R2i)|^) 
and Pl = — P2- Here, R21 is the difference vector between both centres. The propa- 
gation was stopped when the residual Coulomb energy had dropped below 5% of the 
initial value. At the end of each event, the final state was Fourier transformed into the 
momentum space representation. Such a state defines a quantum mechanical (pure) 
ensemble of equally prepared systems. We can extract an arbitrary number of repre- 
sentatives by using the squared amplitude as density function and applying the Monte 
Carlo procedure presented in P] to sample out a number of A^^ correlated pairs. The 2- 
particle correlation function, however, is built by the mixed ensemble which is obtained 
by averaging over the source function. Therefore we have to calculate a large number 
Nf. of events with different initial conditions according to the source densities Eq. (|TUp, 
yielding A'^e ■ Ng pairs which finally obey the correct mixed statistics. 



5 Results 

The source parameters were set to i?s = 5 fm and T = 50 MeV. In order to gain 
an intuitive feeling about the strength of the Coulomb effect, "pions" with different 
charges Z = 1, Z = 3 and Z = 5 were simulated. Since the Coulomb potential 
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between both particles increases quadratically with Z, this set covers a range between 
the reahstic case and the case where the Coulomb energy nearly reaches the same 
magnitude as the localization energy as was discussed in Sect. |1|. = 2000 events 
were generated using a sp- basis set, the average system propagation time was 250 fm/c 
until the Coulomb potential had decayed below 5% of its initial value. Table |I] displays 
some key parameters for the different systems. Clearly the average initial Coulomb 
energy increases quadratically with the particle charge (second column). The third 
column demonstrates that around 70% of this initial Coulomb energy is converted into 
kinetic energy of the wavepackets, independent of the charge. This implies that the 
residual 30% are put into excitation of the p- basis functions, i. e. into deformation of 
the wavepackets. Since with the sp- basis we only have 4 basis functions which are 
orthogonal, we can use the fact that 

\Cs\ + \CpJ + \CpJ + \CpJ = 1 (15) 

to compare the weights that are present in the final state. The fourth column shows 
that in case of Z = 1 in the average only 0.4% of the weight is deposited into the p- 
basis, whereas 99.6% remain in the s- function. This implies that the state remains 
nearly undisturbed. In case of Z = 5 the situation is different: Already 8% of the 
weight is put into the p- functions, in 10% of the events this weight even exceeds 20%. 
A considerable deformation is therefore induced by the Coulomb force. 

For the evaluation of the 2- particle correlation function, each of the 2000 generated 
events was used to extract A^^ = 5000 correlated pairs, yielding a statistics of 10'' pairs. 
Figure |l| displays the resulting correlation functions. Two facts are of importance: 

1. For Z = 1 the distortion of the correlation function is invisible. 

2. For Z = 3 and Z = 5 an increasing deviation from the undistorted correlation 
appears, but the effect is global in the sense that the whole correlation function 
is involved rather than only the part of low momentum differences. 

Several test simulations were performed with different basis sets including sp-, spd- 
and spdf- sets in order to check whether the used sp- basis is already sufficiently large. 
For the Z = 1 simulation no differences were observed, which implies that the sp- 
basis is sufficient to yield reliable results. This is not surprising since the deformation 
of the wavepackets is very small as discussed above. For the Z = 5 simulation some 
dependence on the basis size could be observed. The size of these effects suggests that 
a simulation with spdf- basis instead the sp- basis may lead to further changes of the 
correlation function, which may be as large as 20% of the distortions calculated with 
the smaller basis. Such simulations, however, could not be performed with sufficient 
statistics due to the huge numerical effort. 

As a second time dependent method, the molecular dynamics procedure which was 
discussed in and is based on Bohm's quantum theory of motion, was employed 
to evaluate the pion-pion Coulomb interaction. In this method, only the Gaussian 
s- functions are used. The additional degree of freedom is included by means of an 
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equation of motion for certain test particles, which contains the quantum potential 
and accounts for dispersion- and correlation effects. Since it is the distribution of 
the test particles which is used to deduce the correlation function, the deformation of 
the wavepacket is achieved indirectly due to their Coulomb- distorted positions. This 
ansatz is approximate in the sense that the distorted test particles lead to distorted 
density distributions and therefore to a distorted state, a fact that in an exact treatment 
would change the quantum forces and therefore the particle trajectories and again (in 
higher order) the state. But for only small distortions it was found to be sufficient to 
evaluate the quantum force using the undisturbed state. 

The simulations were repeated for the Z = 1 and Z = 5 case. A number of 
3 ■ 10^ pion pairs were generated for each charge and used to determine the correlation 
functions. For Z = 1 we obtained the same result as with the sp- basis expansion, 
namely no visible Coulomb distortion. Figure ^ only displays the result of the Z = 5 
simulation in comparison to the corresponding result from the basis expansion. An 
overall qualitative agreement is observed, for small momentum differences, however, 
there exist some deviations. This indicates the possible inaccuracy of both methods, 
since for the hypothetical and rather extreme case Z = 5 the Coulomb effects are no 
longer small and both methods begin to suffer from the approximative character of the 
small basis respective the unproper evaluation of the quantum force. For the Z = 1 
simulation, on the other hand, the approximations are well justified and the calculations 
suggest that the pion-pion Coulomb interaction has no visible effect onto the 2- particle 
correlation function and therefore experimental data must not be "corrected" to remove 
a nonexistent effect. 

6 Why scattering theory fails 

Usually, the pion-pion Coulomb effects are treated by considering the time independent 
scattering of plane waves in a Coulomb potential. The shortcomings of this ansatz are 
most obvious if we write the correct quantum mechanical ansatz in the Bohm formalism. 
We may initialize a wavepacket of width cxo with the centre coordinate Ro and a test 
particle at the position Xo and obtain the following equation of motion for the test 
particle: 



with the quantum potential 

Q{t 

which exerts a repulsive force onto the test particle and which when integrated over 
the density distribution yields the zero point energy of the wavepacket |§] . In the plane 
wave limit (cTq oo), the quantum potential vanishes and (|I6|) becomes a classical 
equation of motion. This is the reason why the results from scattering theory can be 



m 



-ViV + Q) 



(16) 



0) 



4mcr? 



(Xq - Rq) 

2a2 



(17) 



6 



made to agree with pure classical trajectory calculations 0. As a crosscheck, we did 
repeat the Z = 1 simulation for wavepackets without dispersion, i. e. for a fixed width. 
The resulting correlation function is in agreement with the predictions of the classical 
"toy model" [^, which proves that it is the dispersion which masks the Coulomb effect. 

In conclusion, we have shown that the correlation functions in pion interferometry 
are not distorted by pion-pion Coulomb effects. Due to the initial conditions, which 
force the pions to be localized within a small volume of a few fm^, the zero point energy 
produces a dispersion which makes the comparably tiny distortion due to Coulomb 
interaction to become unobservable. For heavier particles the strength of the dispersion 
decreases and therefore then the Coulomb effect is expected to be of more relevance. 
We have supported our argumentation by numerical simulations using two different 
methods. In the first (basis expansion), the Coulomb influence is accounted for on the 
level of the state amplitude, in the other (Bohm's molecular dynamics) on the level of 
the equation of motion. Scattering theory fails to produce the correct results, because 
with the removal of the quantum potential in Eq. (|16]) the most important part of 
the dynamical evolution is neglected. This implies that experimental data, which are 
published after "Coulomb correction" , are wrong for small momentum differences. 
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Table 1: Averaged values for the initial Coulomb energy (2'nd column), the difference 
between final and initial kinetic energy (3'rd column) and the degree of excitation of 
the p- basis coefficients (4'th column) for different charges. The numbers (in brackets) 
mean that in 10% of the events that value was exceeded. 



z 


(Eeoui) (MeV) 


(AEkin) (MeV) 




1 


0.21 (0.35) 


0.15 (0.25) 


0.004 (0.01) 


3 


1.9 (3.2) 


1.3 (2.5) 


0.03 (0.09) 


5 


5.2 (9.0) 


3.7 (7.0) 


0.08 (0.2) 
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Figure 1: The 2- particle correlation functions for different charges. The solid curve 
is the theoretical result for noninteracting pions Eq. (|Tl|). 
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Figure 2: Correlation function for Z = 5 and two different numerical methods. The 
solid curve is the theoretical result for noninteracting pions Eq. (pH]). 



